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Abstract 

Employing the formalism introduced by Sinai and Yakhot [PRL, 63(18), p. 1962, 
1989], we study the probability density functions (pdf's) of decaying passive scalars 
in periodic domains under the influence of smooth large scale velocity fields. The 
particular regime we focus on is one where the normalized scalar pdf's attain a 
self-similar profile in finite time, i.e., the so called strange or statistical eigenmode 
regime. In accordance with the work of Sinai and Yakhot, the central regions of the 
pdf's are power laws. But the details of the pdf profiles are dependent on the physical 
parameters in the problem. Interestingly, for small Peclet numbers the pdf's resemble 
stretched or pure exponential functions, whereas in the limit of large Peclet numbers, 
there emerges a universal Gaussian form for the pdf. Numerical simulations are used 
to verify these predictions. 



1. Introduction 

We examine the probability density functions (pdfs) of decaying passive scalars without 
mean gradients under the action of smooth, incompressible and time aperiodic flows in 
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bounded periodic domains. In this situation the evolution of a passive scalar, (l){x,y,t), is 
governed by the advection-diffusion (AD) equation. 

^ + (il • V)0 = «:VV (1) 

Here k represents the molecular diffusivity of the passive scalar and u{x, y, t) is the advecting 
velocity field. The domain (D) under consideration is periodic, specifically we take D to be 
[0, 27r] X [0, 27r] with opposite sides identified (i.e. a 2-torus). 



In this work, we consider time aperiodic velocity fields whose spatial scale of variation 
is comparable to the scale of the domain. Essentially, the flows are of the type encountered 
in chaotic advection^ and the parameters correspond to relatively large Peclet numbers. We 
focus on the case when the normahzed scalar pdf 's attain a self-similar profile. This self- 
similar regime was first described by Pierrehumbert^,^ (who named it a strange or statistical 
eigenmode) and has recently been examined in detail by Fereday and Haynes^ (hereafter FH) 
and also by Sukhatme and Pierrehumbert^ (hereafter SP). The self-similarity of the normal- 
ized pdf is indicative of a non- equilibrium steady state. Physically, as was explained in SP 
and FH, the ingredients in the balance responsible for this state are : (i) a limit on how thin 
filaments can get (as k, > 0) and (ii) the "folding and fiUing" of filaments induced by the 
finite domain. 



To place the things in a proper perspective, we introduce the following scales : L - the 
scale of the domain, /„ - the scale of variation of the velocity field and ls{t) - the maxi- 
mum scale of variation of the scalar field. In terms of these scales the self-similar strange 
eigenmode is characterized by ls{t) ~ ~ L. Due to this similarity of scales the problem 
possesses a global nature (sec FH and SP). Hence, approximations based on scale separation 
(such as shifting to a comoving reference frame), which have yielded excellent results in 
other smooth advection diffusion regimes - sometimes referred to as the Batchelor regime - 
cannot be fruitfully utilized (see^,^ or^ for a recent review). Also, when ls{t) « ~ L (i.e., 

2 



the Batchelor regime), it has been demonstrated that the pdf's are non- universal, i.e. their 
shape evolves in time^. Obviously, in such a situation there is no limiting scalar pdf and the 
theory of Sinai and Yakhot^, which apriori assumes the existence of such a limit, fails. 

Indeed, it is the attainment of a self-similar, i.e. limiting, pdf profile in finite time that 
makes the strange eigenmode regime a suitable candidate for applying the Sinai- Yakhot 
formalism. 



2. The PDF Equation 

Consider the dimensionless normalized variable X = 0/< 0^ >^/^, by assuming the station- 
arity in time of < X^" > for all n, Sinai and Yakhot^ showed, 

(2n - 1) < X^"-^ >=< X^" > (2) 

Qi 

where Qi —< (V0)^ >. Further utihzing Eq. (2), they showed that the pdf of X is given 
by (denoting the sample space variable by the same symbol), 

PiX) = exp[- r-^M (3) 

g{X) Jo g{u) 

where g{X) represents the conditional expectation of the normalized dissipation, i.e 
g{X) —< {'V(f)y/Qi\X >. As it turns out, later work^°,^^,^^,^^ (see^^ for an overview) 
clarified that the pdf of any statistically homogenous twice differentiable random field, say 
ip{x,t), is given by^^,-^^, 

P(V;, t) = , exp[ r (4) 

^ g{ip,t)/Q2 "^70 g{u,t) ^ ^ ' 

Q2 =< {VXf > and r{ii,t) =< (VV)IV' >, gO'^t) =< > represent the condi- 

tional diffusion and conditional dissipation respectively. Furthermore, it was shown that if 
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the moments of il>{x,t) are stationary then-^^/^ ^, 

Ofcourse, Eq. (5) when substituted in Eq. (4) yields a pdf similar to the Sinai- Yakhot 
expression, i.e., Eq. (3). 



A. Conditional Statistics in the Strange Eigenmode Regime 

In the strange eigenmode regime, starting with ls{t = 0) ~ Z„ ~ L, ^ after a transient period, 
it is seen that (see SP and FH), 

< |0(x,y,i)r >~ e-«"* ; t>T (6) 

where T represents the duration of the transient period. More importantly, q;„ = nai, this 
linearity implies the stationarity of the moments of the normalized scalar field ^. Ofcourse, 
given the stationarity, we are justified in using Eq. (5), with X replacing ip. Substituting 
in Eq. (5) from Eq. (6) we get. 



^Ching and Kraichnan^^ hint at the possibility of attaining stationary normalized momemts by 
utilizing a cyclic domain, the self-similar strange eigenmode appears to be precisely this case. 

^Other initial conditions, especially ls{t = 0) « ly ^ L entail an evolution of the scalar field 
through distinct regimes (see SP and FH for details). 

^The exponential decay of moments is also valid when ls{t) « ly, but in this case a„ is a nonlinear 
function of the moment order n, i.e. the moments are not stationary^. 
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Regarding the conditional dissipation, if X and VX are independent, then g{X, t) —< 
(VXy >. This coupled with the fact that g{X,t) is even led Sinai and Yakhot to propose 
the expansion (in the vicinity oi X = 0)^, 

g{X,t) = g{X) =< {VXf > +px' + ...;/? = l^\x=o (8) 
Substituting from Eq. (6), we have (to order X^), 

g(X) = f(l + ?^X^) (9) 

Furthermore, the normalized conditional diffusion and dissipation are, R{X^ — 'r'{X)/ < 
(VXf > and G(X) = g{X)/ < (VX)^ > respectively. Using Eq. (7) and Eq. (9), 

R{X) = -X ; G{X) = 1 + (10) 

Recent numerical work (see FH) suggests that, a2 tends to a non-zero limit as k — > 0, hence 
Eq. (10) gives G{X) — > 1 assuming that /3 does not overwhelm the limit. In other words, 
keeping the assumption regarding (3 in mind, X and VX tend to become independent as 
K — > and we expect the core of the pdf to tend to a universal Gaussian profile. 

Further substituting Eq. (7) and Eq. (9) in Eq. (4) yields (the central part of) the pdf 
of X to be, 

p(x) = C2[i + ^xr^; 7 = i + £^ (11) 

Note that, even though the power-law is in agreement with the work of FH, their arguments 
apply to the tail of the pdf whereas the above expression is valid in the vicinity of X = 0. 
For further elucidation, defining 5 — 2«;/3/q;2, let us examine how the shape of P{X) behaves 
with S. In terms of 6, 

P{X) = C2[l + 5X\^ ; ^ = 1 + ^ (12) 

For large 5 we have 7^1, \a{P{X)) — hi{l-\-5X'^). As both 5 and X are 0(1) quantities, 

all powers of X contribute to ln(P(X)). On the other hand, for small 5 we have 7 1/25 



and ln{P{X)) ~ which is the expected outcome from the earher discussion. Pro- 

files for 5 = 10,0.7,0.001 are shown in Fig. 1. Note that, as 5 decreases P{X) goes from 
resembling a stretched exponential — > pure exponential —>■ (expected) Gaussian function. 



3. Numerical Investigation 

The AD equation was approximated by a lattice map^ followed by diffusion in Fourier space. 
The velocity field is of a single large scale, specifically, we employ the sine flow^^,^^, 

u{x: y, t) = f{t) Ai sm{y + p„) 
v{x, y, t) = {1- f{t)) A2 sm{x + g„) 

(13) 

where f{t) is 1 ior nT < t < {n + l)T/2 and for (n + l)T/2 < t < {n + 1)T. p„, 
(e [0, 27r]) are random numbers selected at the beginning of each iteration, i.e., for each 
period T. Ai.A^ control the strength of the flow. The flow is implemented as a 2D lattice 
map {xn,yn) ^ (^^n+i) Z/n+i)^^)^- A key feature is that the randomness due to Pn,qn breaks 
any barriers which may, and generically do, exist in 2D area preserving mappings^^. 

A typical example : Starting with a mean zero checkerboard initial condition on a 
256 X 256 grid (^i — A2 — 2, k, — 9.33 x 10"^), a typical evolution scenario is shown 
in Fig. 2. As is seen, the normalized moments attain constant values after a transient 
period of about 70 iterations ^. During this transient period P{X) evolves (a remnant of 
the initial double delta pdf can be seen at iteration 10). Finally, after the moments become 

^The transient period may appear large, but it is important to note that the strange eigenmode 
appears only after the scalar field has folded and filled the domain (see SP and FH for details). 

Whereas, significant decay of the variance starts much before this time, specifically when the 
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stationary, P{X) attains a self-similar profile as is seen in lowermost panel of Fig. 2. ^ 

Peclet number dependence : To investigate the effect of changing the Peclet number, we 
run a set of simulations with : (a) fixed fiow strength and varying diffusivity and (b) fixed 
diffusivity and varying fiow strengths. 

• Varying the diffusivity : Keeping Ai,A2 fixed and utilizing the same checkerboard 
initial condition, we vary k. In each case the evolution of the pdf is similar to that 
shown in Fig. 2. The pdfs for k = 2 x 10~^ and k = 5.78 x 10~^ in the self-similar 
stage are displayed in Fig. 3. With respect to the analytical pdf, i.e. Eq. (11), even 
though (3 is an unknown the qualitative similarity between Fig. 3 and Fig. 1 is evident 
(essentially, the dependence of a2 on k is fairly weak, when k changes by an order of 
magnitude, as in the above simulation, a2 changes by a much smaller amount). The 
corresponding plots of the normalized conditional dissipation are shown in Fig. 4. 
Note that for small we have G{X) ~ 1. Also, from Fig. 4 we see that for small k, 
the range of X is quite small, hence in this situation the Gaussian form describes a 
fairly large part of the complete pdf. Next, in Fig. 5 we show the pdfs for a number 
of small diffusivities. Clearly, the core of P{X) tends to a universal Gaussian form. 

• Varying the strength of the fiow : Fixing k = 10~^ we vary Ai,A2. The decay of the 
scalar variance for different fiow strengths can be seen in Fig. 6. Evidently, a2 oc fiow 



diffusive scale is reached. 

^It is important to keep in mind that these results are for large Peclet numbers. In fact, in 
numerical runs with smaller Peclet numbers, even at large times, < X'^^ > fluctuates with a fairly 
large amplitude . 
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strength, therefore for a fixed k, S oc l/(fiow strength). The imphcation being that the 
core of P{X) should tend to a Gaussian function for stronger flows. Fig. 7 shows the 
pdf 's (in the self-similar stage) for two different flow strengths - note the similarity to 
Fig. 3. Furthermore, Fig. 8 shows the pdf 's for a number of simulations with stronger 
flows. Once again, the emergence of a universal Gaussian core is evident. Also, note 
the similarity to Fig. 5. 

4. Conclusion and Discussion 

By applying the formalism introduced by Sinai and Yakhot^ to a decaying passive scalar 
obeying the AD equation in a periodic domain, we obtained an expression for the pdf of 
the normalized scalar field. Broadly categorized as a power-law, the core of the pdf was 
shown to be dependent on the physical parameters in the AD problem. Moreover, we saw 
the emergence of a universal Gaussian core for the pdf in : (a) the limit of small diffusivity 
for a fixed flow strength and (b) the hmit of strong flows for fixed k,. Combining these 
observations we infer the emergence of a universal Gaussian core in the limit of large Peclet 
numbers. Note that, the detailed dependence on Peclet number is not straight forward as 
P{X) is a function of both a2 and k (Eq. (11)) - in turn 0:2 depends on both on k and flow 
strength. Interestingly, for smaller Peclet numbers (i.e. larger 5) the power-law pdf profile 
resembles a pure or stretched exponential function. We believe that this is the reason for 
the mis- identification of pdf profiles in earher work^ (and also SP). 

An intriguing, though poorly understood, feature of the strange eigenmode regime is the 
actual decay rate of the scalar variance - i.e. a2- Note that, in contrast to the present statis- 
tical or strange eigenmode, when the velocity fields are time periodic, and ls{t) 1^ L, the 
scalar field represents a periodic eigenfunction of the AD operator. Hence, the structure as 
well as decay rate of the scalar field are better understood. Regarding these periodic eigen- 
functions, or spatially repeating patterns, see^*^ for experimental results, SP for a physical 
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interpretation,^^ for more recent work,^^ for a mathematically rigorous presentation and^^ 
for an interpretation in terms of the Perron-Probenius operator induced by the underlying 
trajectory problem. Also, see^^,^^ for similar ideas in the case of steady 3D and 2D flows 
respectively. 

Before concluding we would like to put forth a plausible connection between the eigen- 
mode regime and homogenization theory, with the hope of shedding some light on 0:2 ■ 
Broadly, in the realm of homogenization theory it has been possible to show the conver- 
gence, in a coarse grained sense, of the AD equation to a pure diffusion equation for a 
variety of advecting velocity fields (see Section 2 of^^ or^''' for recent reviews). The most 
common situation is when ly « hit), i-©- velocity fields changing rapidly in space but 
combined with either steadiness or periodicity in time^^,^^. The opposite limit, along the 
hues of the work by Kubo^^, is where the velocity fields have longe range spatial correlations 
but change rapidly in time (see section 2.4.1 of^^ or^'''). 

Indeed, it is this second limit where the effective diffusivity of the scalar field is given by 
the Taylor-Kubo formula^^. Recent work has shown that such a diffusive limit exists for a 
broader class of velocity fields, though the formula for the effective diffusivity may not be 
analytically tractable^°,^^. Noting the long range spatial (ly ~ L) and short time (random- 
ness at each iteration) correlations of the velocity field required for the emergence of the 
self-similar eigenmode, we conjecture that the strange eigenmode may be understood as a 
homogenization phenomenon. Hence, averages such as < 0" > obey the diffusion equation, 
in particular a2 can be interpreted as an effective diffusivity akin to the Taylor-Kubo for- 
mula. Not only does this interpretation lend qualitative support to the observed dependence 
of a2 on K as well as the velocity field, it also implies the linearity of with n. 
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Fig. 2. Typical evolution scenario with k = 9.33 x 10~^. Upper panel shows ln(< >) Vs. 
iteration for n=l,2,3 and 4. The middle panel shows the pdf's from iteration 10 to 110 (solid lines 
at iteration 10 and 110) The lowermost panel shows the pdf's from iteration 150 to 350, i.e. when 
the normalized moments have become stationary. 
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Fig. 4. The normalized conditional dissipation {G{X)) Vs. X for the same set of diffusivities as 
in Fig. 3. Note the range of the axes in the two subplots. 
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Fig. 6. Decay of the variance with fixed k and varying flow strengths. Clearly, 02 oc (flow 
strength) . 
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